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I. INTRODUCTION 

Collisions of composite nuclei involve a complicated in- 
terplay and exchange of energy and angular momentum 
between the relative motion and the intrinsic states of 
the nuclei. High precision data for low-energy fusion re- 
actions provide one of the most sensitive tests of such 
interplay. Stationary state coupled-channels descriptions 
have provided a natural methodology to study the ef- 
fects of specific excitation modes of one or both of the 
reactants on the reaction outcomes. The wave function 
in this Schrodinger coupled-channels picture is a linear 
superposition of the states in the model space with a 
definite phase relationship. This coherent linear superpo- 
sition can result in enhancement of the quantum tunnel- 
ing probability and quantum interference. The coupled- 
channels approach has been very successful [l[ in explain- 
ing several collision observables. However, problems re- 
main. Foremost is the inability to describe elastic scatter- 
ing and fusion measurements simultaneously 0, H[ and, 
related, the more recent failure to describe in a physically 
consistent way the below-barricr quantum tunnelling and 
above-barrier fusion yields [H • 

The coupled-channels treatment of collisions of com- 
plex nuclei as closed quantum systems is an approxima- 
tion. In practice, collisions evolve as open quantum sys- 
tems, with innumerable bound and continuum intrinsic 
excitations of the nuclei. In analogy with problems of a 
quantum system in a bath, and following Bohr and Mot- 
telson |5[ , we view the nuclear many-body Hamiltonian as 
a sum of collective, single-particle and coupling terms. In 
a nuclear collision, the collective part comprises the rela- 
tive motion of the nuclei and intrinsic rotational and/or 
vibrational modes. Only a fraction of the Hilbert space 
of this Hamiltonian is used in any feasible full coupled- 



channels calculation. This is because the model space is 
inevitably restricted to selected, most collective intrinsic 
excitations. These collective states define a reduced quan- 
tum system. All other states are weakly coupled to this 
reduced system by residual interactions and constitute an 
"environment" . The key question that arises is: do en- 
vironmental effects influence the reaction dynamics and 
observables, such as angular distributions of products or 
tunneling rates? 

In modelling fusion, the "environment" is assumed to 
come into play only inside the fusion barrier, and is ac- 
counted for in coupled-channels calculations through an 
imaginary potential or an ingoing wave bounday condi- 
tion. However, environmental effects can also be mani- 
fested before the nuclei encounter the fusion barrier, for 
example as real or virtual excitations of giant resonances 
in the individual nuclei by the long-range Coulomb ex- 
citation mechanism. These can be important doorway- 
states to the irreversible loss of kinetic energy (heating 
of the nuclei), as first suggested in Ref. Q for deep- 
inelastic reactions. Another mechanism is complicated 
multi-nucleon transfer channels Q. Measurements have 
shown that deep-inelastic processes occur even at sub- 
barrier incident energies Q , in competition with the pro- 
cess of quantum tunneling, and thus fusion Q. Energy 
loss associated with the deep-inelastic mechanism thus 
could play a significant role in the inhibition of tunneling 
at deep sub-barrier energies. 

The investigation of the effect, on near- and below bar- 
rier fusion, here requires a dynamical model which can 
include both intrinsic state and environmental couplings 
in calculations of the tunneling probability. There are no 
existing realistic theoretical approaches for solving this 
problem. This paper discusses, within a model context, 
(1) how environment-induced irreversibility can be incor- 
porated into the successful coupled-channels framework, 
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and (2) makes a first assessment of its effect on a low- 
energy nuclear collision. The application considered is 
quantum tunneling, relevant to the low-energy nuclear 
fusion hindrance phenomenon. 

The paper is organized as follows. In Section II we give 
a brief survey of theoretical approaches to dissipative dy- 
namics of low-energy nuclear collisions, and discuss the 
suitability of the Lindblad axiomatic theory for the treat- 
ment of energy dissipation on sub-barrier fusion. In Sec- 
tion III we present the coupled-channels density matrix 
approach. Numerical results for our model test case are 
discussed in Section IV. Finally, the summary of the pa- 
per is given in Section V. 



II. THEORETICAL BACKGROUND 

Neither existing models of fusion nor of deep-inelastic 
scattering can address both energy dissipation and quan- 
tum tunneling. The impact of finite lifetimes of excited 
states (e.g., giant resonances) on fusion has been stud- 
ied within a coupled-channels model [1(1 [HI, G3> but 
this approach does not lead to energy dissipation. Di- 
rect damped collisions between complex nuclei have also 
been intensively investigated within various approaches, 
including: (i) transport theories [l3[ based on premas- 
ter, master, Fokker-Planck, Langevin and diffusion equa- 
tions, and (ii) quantum mechanical collective theories 
[Til ]. An appealing semi-classical coupled-channels ap- 
proach combined with a random-matrix model was sug- 
gested by Ko [l5j |. which unifies the statistical and co- 
herent pictures of energy dissipation in deep-inelastic 
collisions. This framework has been succesfully applied 
[i~6l Il7j | to study the excitation of multiphonon giant res- 
onances in heavy-ion collisions at intermediate energy. 
In most of these developments the relative motion of 
the nuclei is described with classical trajectories, whilst 
the coupling to intrinsic degrees of freedom is treated 
either statistically (random-matrix theory) or through 
phenomenological transport coefficients. However, the 
quantum mechanical treatment of the relative motion is 
essential for dealing with quantum tunneling. 

The quantum dynamical model presented in this work 
is based on the time evolution of a reduced density- 
matrix. It provides a consistent description of the transi- 
tion from a pure state to a mixed quantum state during 
the collision. The fundamental equation of motion is the 
Liouville-von Neumann equation for an open quantum 
system, in which a dissipative Liouvillian accounts for 
irreversibility due to interactions of the system with an 
environment. The Lindblad axiomatic approach [Lsl Il9j 
for open quantum systems has been successfully applied 
in nuclear physics, but within rather schematic models. 
For instance, to describe the c harg e equilibration pro- 
cess in deep- inelastic collisions [20j | fission [2l|, decay 
of giant resonances (22j, tunneling through a parabolic 
barrier [23j , and scattering in a two-dimensional inverse 
parabolic potential (24| . These are calculations for a sin- 



gle channel of either one damped oscillator [U [22|, or 
two coupled damped oscillators [2(1 HH . 

In low-energy nuclear collisions, the context of the 
present application, Lindblad's dynamics for the evolu- 
tion of the reduced system is justified (i) because the 
coupling to the complex environment (through excited 
doorway-states and determined by residual interactions) 
is weak, and (ii) because the Markov approximation is 
expected to be valid, the collective motions of the two 
nuclei being slower than the rearrangement of the envi- 
ronmental (nucleonic) degrees of freedom. With increas- 
ing collision energy, to well-above the Coulomb barrier, 
memory effects related to diabatic dynamics [H, [2^, [2?J 
may be important and a non-Markovian Liouvillian may 
be required. The weak coupling between two subspaces 
of the total space of intrinsic nuclear states distinguishes 
between the system and environment. The Lindblad the- 
ory does not require [28| any limitation on the strength 
of the system-environment coupling, although the defini- 
tion of physically well-defined environment states would 
require a careful analysis in the strong-coupling limit. 

An essential effect of the environment on the reaction 
dynamics (unlike the effect of absorptive terms) is to 
progressively destroy the coherent linear superposition 
and the associated phase relationships between different 
channels, introducing quantum decoherence in the sys- 
tem. 

Here, we identify two such (model) sources of deco- 
herence and dissipation. Firstly, an environment inside 
the Coulomb barrier, which is related to the complexity 
of compound nucleus states. Secondly, one with effec- 
tively a long range, associated with decay out of short 
lived (compared to the reaction time) internal vibra- 
tional states, e.g. the giant dipole resonance (GDR) of 
the colliding nuclei, will be shown to be of particular 
importance. The damping of the GDR, due to its ir- 
reversible coupling to a sea of complicated surrounding 
states, which constitute the environment [29(, destroys 
the coherent dynamical coupling with the relative mo- 
tion of the nuclei. Here we show that damping of the 
GDR results in decoherence and energy loss in the re- 
gion where the nuclei overlap, inhibiting tunneling, and 
thus fusion. 



III. COUPLED-CHANNELS DENSITY MATRIX 
APPROACH 

We exploit the time-evolution of a coupled-channels 
density matrix, as is employed in quantum molecular dy- 
namics [30| . The density operator p in Eq. ((T|) is rep- 
resented in an asymptotic (product) basis of states of 
the internal Hamiltonian of the individual nuclei, \i),i — 
1, . . . N (lower indices), and coordinate states describing 
the separation of the two nuclei, \r),r = 1, . . . M (upper 
indices). That is, 

ij,rs 



3 



Crucially, we also add two auxiliary states to the \i) 
basis that allow distinct environmental interactions, as 
described below. The density operator obeys the time- 
dependent Liouville-von Neumann equation 



dp 
dt 



Cp = [C H + Cd]p, 



(2) 



where the total Liouvillian consists of a Hamiltonian part 
ChP = —i[H,p]/h describing the coherent evolution of 
the system with Hamiltonian H, and a dissipative part 
Cd accounting for the interactions with the environment. 
Here, Cd will be assumed to be given by Lindblad's dis- 
sipative Liouvillian [HI, [l^] associated with a Markovian 
semigroup evolution, i.e. 



£dp = ^2(c a pCi - \[cic a ,p} + ) , 



(3) 



where [. . .] + denotes the anti-commutator. Here each C a 
is a Lindblad operator for a dissipative coupling, phys- 
ically motivated according to the specific problem. We 
assume that each coupling a = (Ij) between a given 
state \j) and an environmental state \I) has an associ- 
ated rate T/j, i.e, Cij = [3l| . determined by 
the inverse lifetime of the excited states and the branch- 
ing ratio of its de-excitation, taken to be that when the 
nuclei are well-separated. We note also that: (a) the 
Lindblad Liouvillian has been derived using microscopic 
models [H, H3 , and (b) in contrast to many (dissipative) 
model Liouvillians Eq. ([3]) preserves both the pos- 
itivity and the trace of the density matrix. These are 
essential properties in any realistic application. 

In the model calculations that follow, the basis will 
comprise two asymptotic states (coupled-channels) |1) 
and 1 2) with channel energies ej. Channel |1) is the 
(ground states) entrance channel and is coupled to an 
inelastic state |2) by a coupling interaction Vi2- Two dis- 
tinct sources of irreversibility are also considered, mod- 
elled by two auxiliary (environment) states \X) and \Y). 
The first environmental coupling describes capture by the 
potential pocket inside the fusion barrier. This simu- 
lates the irreversible and dissipative excitations associ- 
ated with the evolution from the two separate nuclei to 
a compound nuclear system. In a stationary states ap- 
proach this loss of flux is approximated by imposing an 
imaginary potential — iW(r), W(r) > 0, or an ingoing- 
wave boundary condition at distances well inside the bar- 
rier. Here, these transitions are described by an auxiliary 
state \X), to which all other states \j) couple, modelled 
by a Lindblad operator C Xj = % /t ff |X)(j|. The ab- 
sorption rate to state \X) is given by 7 rr = W(r)/H where 
W(r) is taken as a Fermi function with depth 10 MeV 
and diffuseness 0.1 fm, located at the pocket radius of 
the nucleus- nucleus potential, ss 7 fm. This choice guar- 
antees complete absorption inside the pocket. The fusion 
probability is defined as the probability accumulating in 
this state \X). 



The second environment, whose explicit treatment will 
be seen to be the most significant at lower energies, is as- 
sociated with the irreversible decay out of intrinsic excita- 
tions of the colliding nuclei. Such decays are independent 
of the dynamical couplings. Specifically, we will associate 
the only excited coupled channel state |2) with the GDR 
excitation. We then introduce a second auxiliary state 
\Y), representing the bath of states in which the GDR is 
embedded, and to which only the GDR excitation \2) is 
coupled. 

Thus, \Y) and/or \X) supplement the two intrinsic 
states |1) and |2) that comprise the two coupled chan- 
nels. Both of the auxiliary states refer to complex ex- 
citation modes of the nuclei, associated with nucleonic 
degrees of freedom and compound nucleus states, respec- 
tively. They provide intuitive and formal channels [33[ 
for describing irreversible coupling and loss of probabil- 
ity from the system to these environments, couplings that 
enter only through the dissipative dynamics term Cd in 
Eq. @. \Y) is also assumed to couple to \X) at the ap- 
propriate range of separations. Probability accumulating 
in state \Y) outside of this \X) pocket may be identi- 
fied with deep-inelastic processes, as will be discussed 
later. These environments and the couplings present in 
the model calculations are represented schematically in 
Figure [1] There, dashed lines indicate regions where the 
intrinsic coupled-channels states |1) and |2) experience ir- 
reversible (environment) couplings to states \X) and/or 

\n 

Upon inserting Eq. |T|) into Eq. @, the following cou- 
pled equations are obtained for the time evolution of the 
density-matrix elements, 



P? 



Explicitly, the Hamiltonian terms are given by 



(4) 



M 



Pu{^ 



A' 



p??(U rr - U ss ) + J2( V ikPkj ~ Pi^J) 



fc=i 



(5) 



where i, j and k run only over intrinsic states |1), |2), 
.... The dissipative terms are given by 



rr rs ' 

kPkk — 2 



k 



prr ^ 
kj/ 



(6) 



where the indices run over all the included intrinsic state 
to auxiliary state couplings and rjj" = +7 rr . In Equa- 
tion ([5]), T, U and V refer to the relative kinetic energy, 
the total bare nucleus- nucleus potential (Coulomb + nu- 
clear) and the coupling interaction between the intrinsic 
channels, respectively. 

While not the technique that is used here, we note that 
this Lindblad dynamical model can also be recast and 
solved numerically within the Monte Carlo wave func- 
tion method, see e.g. Ref. (34J. In that approach, 
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FIG. 1: (Color online) Schematic representations of the dis- 
sipative coupled channels model calculations, with channel 
energies e-,, showing the spatial and energy localization of 
the auxiliary (environment) states |X) and \Y) (shaded rect- 
angles). The upper panel shows the dynamical calculation 
scheme in the presence of environment \X) only. The lower 
panel shows the dynamical calculation scheme in the pres- 
ence of both environments \X) and \Y). Here, dashed lines 
indicate regions where the intrinsic channels |1) and |2) expe- 
rience irreversible couplings to the environmental states \X) 
and/or |Y). 



decoherence and dissipation originate from the intro- 
duction of random quantum jumps in the time evolu- 
tion of the wave function of the system. This unrav- 
eling density-matrix evolution, through stochastic wave- 
function methods [34 ]. shows that the two approaches 
are equivalent and the former takes account of the role of 
fluctuations in the calculation of the expectation values 
and variances of observables. 

The present calculational scheme is based on Eqs. ([5]) 
and © and proceeds as follows. Initially, at time t = 0, 
the nuclei are well-separated, in their ground states, and 
their density matrix describes a pure state with Tr[p] 
= Tr[/3 2 ] = 1. An initial wave-packet describes the rel- 
ative motion of the nuclei. The coupled equations are 
solved numerically using the Faber polynomial expan- 
sion of the time-evolution superoperator [35| . exp(r£), 
and the Fourier method of Ref. [36| for the commutator 
between the kinetic energy and density operator. Having 



solved for the dynamical evolution of the density matrix, 
expectation values of an observable O are now obtained 
from the trace relation (0(t)) — Tr[0/5(£)]. The purity 
of the initial state, conserved under Hamiltonian unitary 
evolution, will be destroyed (Tr[/5 2 ] < 1) if the environ- 
ment causes a loss of quantum coherence. This decoher- 
ence can thus be quantified via this loss of density-matrix 
purity, or equivalently by an increase of the linear entropy 
1-Tr[p 2 ]. 



IV. NUMERICAL RESULTS AND DISCUSSION 

So as to make contact with the coherent, Schrodinger 
picture, the model Hamiltonian we used was chosen to 
coincide with that of the coupled-channels fusion model 
CCFULL [33] • Specifically, our model calculations use 
physical parameters relevant to the 16 + 144 Sm reac- 
tion at collision energies below its nominal fusion bar- 
rier, Vb = 61.1 MeV. We assume zero relative orbital 
angular momentum between the reactants. The form of 
the bare nuclear potential between the two nuclei, con- 
sistent with the stated Vb, is a Woods-Saxon potential 
with (V ,r ,a ) = (-105.1 MeV, 1.1 fm, 0.75 fm). The 
Coulomb potential was that for two point charges. The 
16 projectile was taken to be inert and the 144 Sm tar- 
get was allowed to be excited to a GDR vibrational state. 
The dynamical nuclear coupling of the ground state 1 1) to 
the vibrational state |2), with excitation energy E x - = 15 
MeV, has a macroscopic deformed Woods-Saxon form 
with a deformation parameter of j3\ = 0.2. 

The time step for the density-matrix propagation was 
At = 10~ 22 s, and the radial grid (r — — 250 fm) was 
evenly spaced with M — 512 points. The relative mo- 
tion of the two nuclei in the entrance channel |1) was de- 
scribed by a minimal- uncertainty Gaussian wave-packet, 
with width 00 = 20 fm, initially centered at r = 150 fm, 
and was boosted towards the target with the appropriate 
average kinetic energy for the entrance channel energy Eq 
required. The FWHM energy spread of the wave-packet 
is ~ 3%. The numerical accuracy of the time evolution 
was checked using a fully coherent, time-dependent cal- 
culation, excluding coupling to states \X) and \Y). It 
was confirmed that the normalisation and purity of the 
density-matrix, Tr[/5] — Tr[p 2 ] = 1, and the expectation 
value of the system energy Tr[Hp] were maintained with 
high accuracy over the required number of time steps, 
typically 700 for the full duration of the collision. 

The importance of the two, spatially distinct, sources 
of environment couplings were studied. Calculations were 
first performed in the scheme shown in the left panel of 
Figure [1] Here the intrinsic coupled channels |1) and |2) 
also couple to the capture state \X). Calculations were 
carried out for E = 45, 50, 55 and 60 MeV incident en- 
ergy. The calculated state purity Tr[/5 2 ] and the energy 
dissipation Tr[H (po — p)] post the collision (after 700 time 
steps) are shown in the left panels in Table HI For suffi- 
ciently sub-barrier energies, Eq < 55 MeV, it is evident 
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TABLE I: The calculated density matrix purity Tr[p 2 ] and 
energy loss AE =Tt[H(po — p)] following time-evolution (for 
700 time steps) when including only the state | X ) (left entries) 
and both states \X) and \Y) (right entries) environmental 
couplings. The GDR coupling strength used was f3\ = 0.2. 





State \X) 


States \X) and \Y) 


E (MeV) 


Tr[p 2 ] 


AE (MeV) 


Tr[p 2 ] 


AE (MeV) 


45 


1.0000 


0.0004 


0.9196 


1.8718 


50 


1.0000 


0.0004 


0.8977 


2.6744 


55 


0.9996 


0.0109 


0.8759 


3.6100 


60 


0.6067 


14.862 


0.5127 


18.908 



that time-evolution in the presence of state \X) essen- 
tially maintains coherence and is non-dissipative. There 
is however loss of purity and dissipation at the highest 
energy. It is interesting therefore to compare the density- 
matrix and the Schrodinger predictions of CCFULL (that 
uses an ingoing wave boundary condition). This is done 
here only for calculations of the tunneling probability 
P(Eq), in a relative s-wave, shown in Figure^a). These 
comparisons, of necessity, require convolution of the £ = 
partial wave penetrabilities T (E) from CCFULL with the 
energy distribution f(E, E ) of the chosen initial wave 
packet. That is, P(E ) = J dE f {E , E ) T (E) . The 
P(Eq), shown as a function of Eq/Vb in Figure Ela), 
are in very good agreement showing the appropriateness 
of stationary state coupled-channels calculations for this 
observable within the dynamical scheme of states |1), |2) 
and \X). It is our contention that the dissipation as- 
sociated with state \X), while significant at 60 MeV, is 
strongly localised inside the barrier and thus does not 
impact upon the barrier penetrability. We will now show 
that the same is not true for the more spatially-extended 
dissipation due to the GDR decay environment \Y). 

The treatment of the irreversible GDR decay (with 
a spreading width of 6 MeV) to the bath of surround- 
ing complex states (represented by \Y}) was included by 
switching on the coupling of the intrinsic inelastic state 
1 2) to \ Y) . This is the dynamical scheme of the right hand 
panel in Figure [TJ Unlike the coupling to \X), a major 
part of the inelastic excitation of the system gives access 
to \Y) before the wave packet encounters the fusion bar- 
rier. The onset of decoherence, the purity of the density 
matrix, and the associated energy dissipation are shown 
in the right hand entries in Table [J and, as a function 
of time evolution, in Figure 2(b), the latter for Eq = 45 
MeV. 

Summing over all internal states of the density ma- 
trix gives the total diagonal elements in coordinate space, 
which represent the wave-packet at a given time. Snap- 
shots of the wave-packet in the interaction region (for Eq 
= 45 MeV) are shown in Fig. 3. The curves are shown for 
times t = (dashed; the initial state), 400 x 10~ 22 s (dot- 
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FIG. 2: (Color online) (a) The energy dependence of the s- 
wave tunneling probability calculated with the density ma- 
trix (solid points) and the coupled-channels CCFULL methods 
(full line), (b) The time evolution of the density matrix pu- 
rity Tr[/3 2 ] and the energy loss Tr[_ff(po — p)\ with decoherent 
\X) and \Y) state dynamics: corresponding to the right-hand 
E = 45 MeV entry in Table 1. 



ted; near to the time of closest approach) and 700 x 10~ 22 
s (thin solid; post the collision). The Figure shows the 
results from (a) the coupling to state \X) , and (b) to 
both \X) and \Y). When the wave-packet tunnels into 
the pocket (from dashed to dotted lines), the short range 
coupling to \X) leads to trapping of flux from |1) and 
1 2) inside the potential pocket. This reveals itself as 
an unchanging probability for radii r < 7.5 fm as the 
main body of the wave-packet leaves the interaction re- 
gion (from the dotted to the thin solid lines) in Figure 
3(a). The additional effect of turning on the coupling 
between states |2) and \Y) is to trap probability under 
the barrier, as is shown in Figure 3(b). This reduces the 
component of the wave-packet that reaches the potential 
pocket, inhibiting the quantum tunneling. 

The probability trapped under the fusion barrier is as- 
sociated with GDR collective vibrational energy being ir- 
reversibly removed from the coherent dynamics into innu- 
merable channels (heat) . This is then no longer available 
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FIG. 3: The inter-nuclear potential (thick curve) and the time 
evolution of the wave-packet for ie O + 144 Sm at Eq = 45 
MeV: (a) including only coupling of auxiliary state \X), (b) 
including coupling of the states \X) and \Y). The wave packet 
is plotted at times t = (dashed curve), 400 x 10" 22 s (dotted 
curve) and 700 x 10~ 22 s (thin solid curve). 



for relative motion, or tunneling. Such energy loss can 
be correlated with deep inelastic processes, seen experi- 
mentally, that compete with fusion in reactions involving 
heavy nuclei 

Fig. fjja) shows the time evolution of the probability 
trapped in the potential pocket, state \X), for E = 45 
MeV. We comment that, when including the inelastic 
channel \2) but not \Y), the nucleus-nucleus potential 
renormalization leads to the expected enhanced penetra- 
bility from the inelastic channel coupling, compared to 
the purely elastic (|1) plus \X)) calculation. The deco- 
herent dynamics due only to environment \X) gives the 
(full curve). By comparison, the calculation that also 
includes the GDR doorway-state decay to |V) leads to 
a suppression (dotted curve and arrow) of the popula- 
tion of state \X). Additional irreversible processes other 
than excitation of the GDR are also likely to contribute 
to the deep-inelastic yield, such as complicated multi- 
nucleon transfers Q- To simulate these very simply, the 
assumed state |1) to \2) coupling strength was increased. 




0.2 0.4 0.6 0.8 1 



Pi 

FIG. 4: (Color online) (a) Time-dependence of the probability 
trapped in \X) for Eo — 45 MeV. The full curve includes 
states |1), |2) and \X). The dotted curve adds the irreversible 
decay of |2) to \Y). The calculations are for (5\ = 0.2. (b) 
Calculated suppression of the probability trapped in \X) as a 
function of the assumed /3i value for Eo = 45 MeV (dashed 
curve) and 55 MeV (solid curve). 



Fig. 0Jb) shows the dependence of the calculated tunnel- 
ing suppression on the assumed (3\ strength for Eq = 45 
(dashed curve) and 55 MeV (solid curve), where we note 
that larger j3\ result in both an increase in the strength 
and the range of the coupling formfactor to the inelastic 
state 1 2). 

V. SUMMARY 

A quantum dynamical model based on time- 
propagation of a coupled- channels density matrix has 
been presented and is shown to describe the transition 
from pure state (coherent) to mixed state (decoherent 
and dissipative) dynamics during a nuclear collision. The 
calculations exhibit both decoherence and energy dissi- 
pation and so go beyond coherent coupled-channels ap- 
proaches. Decoherence, originating here from the irre- 
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versible decay of a giant-dipole vibrational state of the 
heavy target nucleus to surrounding states, is shown to 
result in hindrance of quantum tunneling. Developments 
of the model calculations to include (a) non-zero relative 
orbital angular momenta between the reactants, (b) ad- 
ditional intrinsic channels, and (c) a more detailed con- 
sideration of other processes, such as multi-nucleon or 
cluster transfer reactions, are necessary in order to con- 
front experimental measurements. 
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